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Abstract 

At large parallel magnetic field Bu, the ground state of bilayer quantum 
Hall system forms uniform soliton lattice phase. The soliton lattice will melt 
due to the proliferation of unbound dislocations at certain finite temperature 
leading to the Kosterlitz-Thouless (KT) melting. We calculate the KT phase 
boundary by numerically solving the newly developed set of Bethe ansatz 
equations, which fully take into account the thermal fluctuations of soliton 
walls. We predict that within certain ranges of Bn , the soliton lattice will melt 
at Tkt- Interestingly enough, as temperature decreases, it melts at certain 
temperature lower than Tkt exhibiting the reentrant behaviour of the soliton 
liquid phase. 
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When the interlayer spacing d is comparable to the mean particle distance, the bilayer 
quantum Hall system at total filling factor u tot = 1 can exhibit the quantum Hall effect due 
to the strong interlayer correlations [0-0. In the presence of large parallel magnetic field 
B\\, the ground state of the system is known to form a uniform soliton lattice (SL) phase 
made of the periodic array of the phase solitons of the field variable 9(r), which represents 
the relative phase difference between electrons in two layers ||. Using the isospin language 
§|, the effective energy functional for the bilayer system can be written by 

E{6)= \ Sr {^ Ps|W|2 "2^ cos ^( r ) + ^]}> (!) 

where p s is the isospin stiffness, t = toe'® 2 ^^ the interlayer tunneling amplitude, £ = 
(hc/lelB^) 1 / 2 , and Q = 2ir/L\\ with Ln = &o/Bnd, which defines the length associated with 
one flux quantum $ enclosed between two layers. We set the magnetic length 1=1. This 
continuum model is valid, when the field 8(r) varies smoothly over the lengths £, L\\ [|J. With 
the increase of B», the system exhibits a quantum phase transition from commensurate (C) 
to incommensurate (I) phase at Q c = {2/ r K){2t/ r Kp s ) 1 / 2 at T = @,[|. The incommensurate 
phase at Q > Q c describes the uniform SL phase at zero temperature. 

In the paper, we study the thermodynamics of the SL phase through a mapping of the 2D 
statistical mechanical model in Eq. (|l|) to the ID quantum sine-Gordon (QSG) Hamiltonian. 
When one neglects the compactness of the field variable 9, which will be reasonably valid for 
T < (7r/2)p s , the ID QSG model is exactly soluble by the Bethe ansatz method. The Bethe 
ansatz solution provides the thermodynamic CI phase boundary Q C (T). Depending on the 
values of C = t/(32np s ), the critical value Q C (T) will increase (decrease) with temperature 
for C > 1 (C < 1). Within the incommensurate phase, the ground state of the system is the 
uniform SL phase at T = 0. As T increases to about (n/2)p s , the compactness of the angle 
variable 9 will become important, which introduces dislocations to the system leading to 
the KT melting of the SL phase. Recently Hanna et al. have studied the melting of the SL 
phase by analyzing the elastic moduli based on the rigidity of the zero temperature ground 
state However, the renormalizaion of the elastic moduli due to thermal fluctuations 
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of soliton walls has not been taken into account. Furthermore, as one approaches to the CI 
phase boundary, the soliton density will vanish. In order to understand the melting of the 
SL phase near the CI phase boundary, it is crucial to include the thermal fluctuations and 
the soliton density variations as well. 

In the paper, we develop a new set of Bethe ansatz equations, which fully take into 
account both the thermal fluctuations of the soliton walls and the density variations. By 
numerically solving the Bethe ansatz equations, we have calculated both the compression 
modulus k xx and the shear modulus n yy of the SL phase. By making an asymptotic expansion 
near the CI phase boundary, we have shown that k xx goes like (2-r/M) 1 / 2 (Q — Q C (T)) 1 ^ 2 p s 
and K yy = (8Mr 3 ) 1/2 (Q - Q C (T))~ 1/2 p s , where M is the soliton mass which varies with 
t and r = T/(8np s ). This asymptotic behaviour explicitly confirms the predictions made 
by Coppersmith et al. J7j. Based on the elastic moduli thus obtained, Tkt is calculated 
and compared with the zero temperature estimates |]|§. We predict the following reentrant 
behaviour of the soliton liquid phase. At certain values of Q slightly below Q c (0), the system 
initially stays at the C phase. As T increases, it makes a transition to the soliton liquid 
phase. With the further increase of T, the soliton liquid solidifies to the SL phase due to the 
rapid increase of the soliton density, which subsequently melts reentering the soliton liquid 
phase. 

Following the effective energy functional of Eq. ([[]), the low temperature thermodynamics 
of the system can be defined by the statistical partition function Z = f[D9]e~ E ( 9 ^ T . One 
can map the statistical ensemble summation into quantum transfer matrix by identifying 
Z = lim^oo Tr 



e -RH 



H = dx 



W? + %r (dJ? -^f cos{e + Qx) 



(2) 



If one neglects the compactness of the angle variable #(r), the canonical conjugate mo- 
mentum II(r) can be easily integrated out leading to the corresponding classical partition 



function. We make the following change of variables 0(r) = y p s /T [6(r) + Qx], which leads 
to the quantum Hamiltonian of the ID sine-Gordon model with external U(l) coupling, 



H 



dx 



\li 2 + Ud x <pf - A^d x cf> + 2/icos {{3<f>) + ^-p s Q 2 

L L 111 21 



(3) 



Here fi = t/(4irT), [3 = / p s , and A = 2nQp s /T. The compactness of 6>(r) would 
restrict Il(r) to the integer values. Hence the ground state energy £ of the ID QSG model 
corresponds to T '/T ', where T is the free energy of the 2D statistical mechanical system. 

The ground state of the ID QSG model is exactly soluble by the Bethe ansatz method. 
There is a competition between the finite soliton mass M which prefers the commensurate 
phase and the external field A coupled to the topological charge which favors the incom- 
mensurate phase. This is a 'quantum' version of the CI transition [[|, where the soliton 
mass includes full 'quantum' fluctuation effects. With external U(l) field, a soliton can be 
created with energy cost AS = M cosh — A. The soliton mass M is given by 



M 



2r(§) 



7T/i 



(4) 



where p = T/(8np s — T) and represents the rapidity of the soliton [|T0[] . 

For A < M, the ground state will be the vacuum: the commensurate phase. However 
for A > M, a nontrivial vacuum will arise, since the energy cost to create a soliton can be 
negative. In this case, the ground state will be described by the soliton condensations: The 
SL phase. The CI transition will occur at A = M. By equating the two quantities, one can 
obtain the exact CI phase boundary Q C {T) 



Qc(T) 



c- 



Tfl + r) 



2(1-t) 



(5) 



where we have neglected the Q-dependence of t, since the continuum model is only valid 
for Q < 1. The parameter C is given by t/(4To). The reduced temperature variable r is 
between and 1. As r — > 0, one can reproduce the classical limit: Q c (0) = IQ\/C/it. 

Now we want to study the melting of the SL phase by numerically solving the Bethe 
ansatz equations. Using the fact that the exact scattering matrix is known for the ID QSG 
model, one can calculate the ground state energy and the soliton density n s (Q,T) as a 



function of Q and T [10 . The Bethe ansatz equation is given by 



27rp(6) = Mcoshe + [ B d&(p(Q - e')p(e'), (6) 

where p(0) is the density of solitons between and + dO per unit length and the integral 
kernel ip(Q) is given by 



oo 

<p(Q)= I dwe^ e ? \ A - ' v . (7) 

J-<*> 2cosh(f^J sinh(^) 

In terms of the density function, one can calculate the ground state energy as follows 

AS = [ B d0(Mcosh0 - A)p{G), (8) 

where the ground state is determined by imposing the condition: dAS/dB = 0. By making 
an asymptotic expansion of Eq. (^) near the CI phase boundary, one can confirm that the 
soliton density exhibits a power- law behaviour: n s = (2M/t) 1 ' 2 (Q - Q C (T)) 1/2 /(2tt) §. 

As T increases to about (n/2)p s , the compactness of the variable 8(r) becomes important, 
which introduces topological defects (dislocations) into the system leading to the KT melting 
of the SL phase at temperatures much below the CI transition temperature. Although the 
ID QSG Hamiltonian does not involve dislocation excitations, one can get a reasonable 
estimate for Tkt by analyzing the elastic moduli of the SL at finite temperature. The elastic 



sinh ^p-V"' 
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2 cosh i 




! sinh ( 


Trpuj \ 
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moduli Kji j £1X6 defined as follows 



d 



2 



K 



lJ dQidQi 

where Q x s is given by 2n/L s , where L s is the lattice constant of the SL phase along the x 
direction, Q y s = tan Q x s with the tilt angle of the soliton walls, and i, j = x, y. In order to 
obtain the elastic moduli of the system, one needs to calculate the energy of the system by 
varying the soliton density from the ground state value and/or tilting the solitons from the 
vertical orientations. We first calculate the compression modulus k xx . By shifting the fermi 
momentum from B to B + e, one can increase (decrease) the soliton density for positive 
(negative) values of e. The soliton density p(0) is a function of both e and 0. For small e, 
one can expand the soliton density up to quadratic order in e 



■F, (9) 



e 2 



p(e)^p(e) + ep l b (e) + - P 'l(e), (10) 

where p(0) is the soliton density profile of the ground state, p' b = dp/de\ e=0 , and p' b ' = 
d 2 p/de 2 \ e= Q. At the lowest order of e, we obtain the usual Bethe ansatz equation for the 
ground state as shown in Eq.(^) with the condition dAS/de = 0. Up to the order of e 2 , we 
obtain the following two equations 

2^(6) = b(e - b) + v {q + B)}p(B) + [ B dev(0 - 0'K(©'), (11) 

J — B 

dp(B) 



2 V b '(0) = ^ [<p(Q + B) - <p(Q - B)} p(B) + [^(0 - B) + <p(B + B)] 



dB 

+ f d0V(0-0'K(e')- (12) 

J —B 



The energy of the system is given as follows: AS = AS + \Xb, where AS is the ground 
state energy and Xb is given by 

fB 



Xb = 2M sinh Bp{B) + 2(M cosh B — A) 



-QB- + 2p 6 (5) 



+ / d0(Mcosh0- A)pl{&). 

J — B 

(13) 



The Q x can be calculated as follows: Q x = 2ic J?gl e dQp(Q) = Q x + a b e + 0{e 2 ), where 
Q x = 2tt J B B dQp(&) is the ground state value and a b = 4irp(B) + 2tc J^ B dQp' b (&). The 
compression modulus k xx is given as follows: d 2 T/dQ x2 = Txb/al 

Now we calculate the shear modulus K yy . To obtain the shear modulus, one needs to 
fix Q x and vary Q v s alone. Global shift of the rapidity will not only tilt the soliton walls 
but also change the soliton density along the ^-direction. In order to vary Q y s alone, one 
needs to change the fermi momentum B to B + e and rapidity by rj simultaneously along 
the following trajectory: e = — a^rj 2 for some constant ax ||llj| - By making a perturbative 
expansion in 77, we obtain the following two Bethe ansatz equations up to r] 2 , 

2V S (0) = b(© - B) - <^(0 + B)]p(B) + f B deV(6 - &)£{&), (14) 

J — B 
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2 V;(©) = qq Me + b)- <p(e - b)} p(B) + (<p(e -b) + <p(e + b)) 

+ / B deV(e-e'K(e'), (is) 

where one can easily notice that p' s (Q) is an odd function of 9 and p"(0) an even function. 
The energy of the system along the trajectory is given by 

AS = AS + ^Xs, (16) 

where \s is given by 

+ ( B rf6(Mcoshe- A)p"AQ). 

J-B 

(17) 

Since Q y s = tan0Qf = vQsi the shear modulus K yy is given by the following rela- 
tion: d 2 TjdQ y2 = Txs/(Qs) 2 - The off-diagonal elements of «jj are zero. Using the 
asymptotic expansion, one can show that near the CI phase boundary, k xx goes like 
(2r/M) 1 /2 (Q _ Q c (T)) 1/2 p S) and K yy = (SMr 3 ) 1 / 2 (Q - Q c (T)y 1/2 p s |TJ. 

In Murphy et a/.'s experiment, for equally populated layers, the experimental values 
of the various parameters are given as follows: p s = 0.35K,to = 1.2K,£ = 126 A, and 
d = 200 A [Q. For the above sample, the parameter C can be estimated to be about 
0.033. In the inset of Fig. (1), the open circles represent k xx and the open squares n yy 
at fixed value of Q = 0.919, which obey the correct asymptotic behaviour near the CI 
phase boundary. We have confirmed that both k xx and K yy approach to p s at large Q 
fPlj| . The KT transition temperature Tkt can be estimated by solving the following equa- 
tion: k B T/(%/2)p s = [K xx (T)K, yy (T)] l l 2 / p s 0. In Fig. (1), the closed circles represent 
(k>xx(T) K yy (T)Y I 2 1 '(16 p s ) as a function of r. The intersections with the solid line with slope 
1 locate the positions of Tkt- It can be shown that for Q m i n < Q < Q c (0), the equations 
have two solutions, where Q m i n = 0.9, Q c (0) = 0.927. In Fig. (2), we plot the KT melt- 
ing temperature as a function of Q. The closed circles represent the KT phase boundary 
from our Bethe ansatz calculation, the open squares from Hanna et al. based on the zero 
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dp(B) 
dB 



Xs = 2M sinh Bp{B) + 2(M cosh B — A) 



dp(B) 
dB 



temperature value of the elastic moduli, and the solid curve the CI phase boundary. For 
Qmin < Q < <5c(0), the system initially stays at the C phase. As T increases, it makes 
a transition to the soliton liquid phase. With the further increase of T, the soliton liquid 
becomes the SL phase due to the rapid increase of the soliton density upon entering the 
CI phase boundary. Subsequently it melts exhibiting the reentrant behaviour of the soliton 
liquid phase. We have confirmed that at large Q, Tkt approaches to (n/2)p s [fLl] . Read 
argues that at T = 0, quantum fluctuations are not important, since the domain 'sheets' are 
marginally rough P,12|. However at finite T, thermal length scale £p = {hv/ksT) becomes 
finite. Hence near the CI phase boundary and/ or KT phase boundary, the large distance 
thermal fluctuations will be still important. 

To summarize, we have studied the melting of the SL phase by numerically solving the 
newly developed set of Bethe ansatz equations, which fully take into account both the ther- 
mal fluctuations of the soliton walls and the density variations. Based on the elastic moduli 
thus obtained, Tkt is calculated. We predict that the system will exhibit the reentrant 
behaviour of the soliton liquid phase within certain ranges of Q. Our proposed KT phase 
boundary can be experimentally observed by measuring the longitudinal resistance of the 
high purity bilayer sample. Towards the completion of the paper, we became aware that 
the preprint cond-mat/ '0201343 also studied the bilayer quantum Hall system based on the 
exact solution of the ID sine-Gordon model. In our paper, we mainly focus on the role of 
the compactness of the angle variables and explicitly calculate the KT phase boundary. 
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FIGURES 

FIG. 1. Determination of the KT melting temperature. The closed circles represent 
(Kxx(T)Kyy(T)) 1 / 2 1 (16/9 S ) as a function of r. The intersections with the solid line with slope 1 
locate the positions of Tkt- In the inset, the open circles represent the compression modulus n xx 
as a function of r at fixed Q = 0.919 and the open squares the shear modulus K yy . 

FIG. 2. The KT melting temperature as a function of Q. The KT phase boundary is shown 
and compared with the result of Hanna et o/.[6], where tkt = 16r. The closed circles are from 
our Bethe ansatz result and the open squares from Hanna et ai. The solid line represents the CI 
phase boundary. 
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